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AO Vel: The role of multiplicity in the development of 
chemical peculiarities in late B-type stars* 
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ABSTRACT 

We present high-resolution, high signal-to-noise UVES spectra of AO Vel, a quadruple 
system containing an eclipsing BpSi star. From these observations we reconstruct 
the spectra of the individual components and perform an abundance analysis of all 
four stellar members. We found that all components are chemically peculiar with 
different abundances patters. In particular, the two less massive stars show typical 
characteristics of HgMn stars. The two most massive stars in the system show variable 
line profiles indicating the presence of chemical spots. Given the youth of the system 
and the notable chemical peculiarities of their components, this system could give 
important insights in the origin of chemical anomalies. 

Key words: stars: binaries - stars: chemically peculiar - starsdndividual: HD 68826 



1 INTRODUCTION 

AO Vel (=HD 68826) was classified as BpSi by 
iBidelman fe MacConneU Jl973l) and was reported as 
an ec lipsing binary by IClausen. Gimenez. fc van Houtenl 
In fact it is one of the only two double-lined 
eclipsing binaries with a BpSi component known to date 
l|Hensberge et al.ll2004). Fr om li ght-time effect on the times 
of minima, Clausen et al. I (|l995l) deduced the presence of a 
third body. 

In our previous paper (|Gonzalez et al.l I2006T ). using 
FEROS spectroscopic observations, we discovered that this 
system is actually a spectroscopic quadruple system with 
components close to the ZAMS. The four stars form two 
close spectroscopic pairs (periods of 1.58 and 4.15 days) 
bound gravitationally to each other in a wide eccentric orbit 
with a period of 41 yr. In that paper we combined our radial 
velocity (RV) measurements with the available photomet- 
ric data to derive orbital parameters for both binary systems 
and to calculate the absolute parameters of the eclipsing sys- 
tem. For the first time, direct determination of the radius 
and the mass was obtained for a BpSi star. 

In this work we present high-resolution, high signal-to- 
noise UVES spectra, which were used to perform an abun- 
dance analysis of all four components of this multiple sys- 
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tem. In Sec. 2 we present the observations and describe the 
reconstruction of the component spectra. In Sec. 3 we an- 
alyze the spectral characteristics of each component and 
present the results of the abundance analysis. In the last 
Section we discuss the main results and the occurrence of 
chemical peculiarities in binary and multiple systems. 



2 OBSERVATIONS AND SPECTRA 
RECONSTRUCTION 

Three spectra were obtained in service mode with UVES 
at VLT-UT2 telescope in October 2005. The spectra have 
been taken with the 0.4arcsec slit for the blue arm and the 
0.3arcsec slit for the red arm on three consecutive nights 
with two different dichroics to achieve the highest UVES 
resolution of 110,000 in the red spectral region and 80,000 
in the blue spectral region. We used exposure times of 20-30 
min, obtaining a S/N ratio above 200 in the spectral range 
4000-8000A. 

These spectra are analyzed here along w ith five FEROS 
spectr a described in the previous paper by IGonzalez et al.l 
(2006). To calculate separate spectra for the four compo- 
nents of the system and t o measure their RVs, the iter- 
ative method described by IGonzalez fc Levatol (2006) was 
adapted for the multiple system AO Vel. This algorithm 
computes the spectra of the individual components and the 
RVs iteratively. In each step the computed spectra are used 
to remove the spectral features of all but one component 
from the observed spectra. The resulting single-lined spec- 
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Table 1. Modified Julian dates of UVES observations, orbital phases, measured RVs, and signal-to-noise ratio at three different wave- 
length. 



MJD 


Phase AB 


Phase CD 


RV A 




RV C 




RV D 




S/N A42 oo 


S/Na5600 


S/Na8500 








krns" 1 


km s _1 


km s — 


i 


krns" 


1 








53669.3357 


0.7133 


0.8722 


172 ± 4 


-180 ± 4 


88.2 ± 


2.2 


-49.7 ± 


0.9 


250 


300 


175 


53670.3346 


0.3436 


0.1129 


-149 ± 3 


172 ± 5 


-38.6 ± 


2.2 


88.2 ± 


0.9 


310 


375 


190 


53671.3523 


0.9859 


0.3582 


52 ± 8 


-9 ± 8 


-50.2 ± 


4.7 


102.2 ± 


1.6 


235 


320 


180 



tra are used to measure the RV of that component and to 
compute its spectrum by combining them appropriately. 

In these calculations we combined the three new UVES 
spectra with our four FEROS observations taken out of 
eclipse. Since the resolution and S/N-ratio of the UVES 
spectra were much higher than for FEROS spectra, we as- 
signed higher weight to these spectra in the last iterations. 

Table 1 lists the measured RVs along with the phases 
(zero at conjunction) corresponding to the orbits AB and 
CD. For the calculation of orbital phases we have applied 
a time correction of +0.0236 d and -0.0337 d for the system 
AB and CD respectively, in order to take into account the 
light time effect. The last spectrum in this table was taken 
at a phase of partial eclipse of the primary star A. We have 
taken into account this fact during the calculations of the 
component spectra, considering that the light contribution 
of A is smaller in that spectrum and therefore the line in- 
tensities are expected to differ from those in the remaining 
spectra. If, for out-of-eclipse phases, star A contributes a 
fraction l a to the total light of the system, and during a 
phase of partial eclipse the light of this star is reduced by a 
factor /, then the apparent intensity of the spectral lines of 
the star A is reduced by a factor f /(f -la + ^ — la), and that of 
the remaining components is increased by • l a + 1 — la)- 
These scaling factors were estimated to be, in the case of 
our third UVES spectrum, 0.85 and 1.10, respectively. 

We note that the amplitude of the RV curves of the sys- 
tem AB might be affected by the asymmetry of the spectral 
line profiles, which present variations for several elements 
(see next Section). This error source has not been considered 
in the formal error quoted in Table 1 , since it is difficult to 
be estimated without determining the surface distribution 
of the various atomic species giving rise for spectral lines 
used for RV measurements. 

The obtained new RV measurements were in general 
agre ement with the orbit s published in our previous pa- 
per ([Gonzalez et al.ll2006l l. especially for the eclipsing bi- 
nary AB, for which some orbital parameters are fixed by the 
photometric data. In the case of the pair CD, the new data 
allowed a significant improvement of its orbital parameters. 
The orbital parameters recalculated combining all available 
RV measurements are listed in Table 2. In Fig. 1 we present 
the RV curves for the eclipsing binary system AB and for 
the less massive system CD. In the case of the binary AB, 
which presents apsidal motion, the solid line corresponds to 
the orbit calculated for the epoch of FEROS observations 
and the dotted line to the epoch of UVES observations. 

Once the spectra of the stellar components have been 
reconstructed, they must be scaled to recover the intrinsic 
intensities of the spectral lines in those stars. To this aim, 
as in our previous paper, we used the photometric relative 



Table 2. Orbital parameters for the spectroscopic binaries AB 
and CD. 



Parameter 


Unit 


Value 




P C D 


days 


4.14933 


± 


0.00004 


T CD (conj) 


MJD 


53300.539 


± 


0.007 


V~,CD 


km s _1 


21.5 


± 


0.4 


K C 


km s" 1 


94.6 


± 


1.2 


K D 


km s~ x 


103.9 


± 


0.7 


CCD 


Rq 


0.017 


± 


0.006 


<+>CD 


deg 


261 


± 


22 


acD sini C D 


Rq 


16.48 


± 


0.12 


QCD 




0.917 


± 


0.016 


M c sin 3 i C D 


M e 


1.82 


± 


0.03 


M D sin 3 icd 


M@ 


1.67 


± 


0.05 


V 1 ab 


km s _1 


8.1 


± 


1.8 


K A 


km s _1 


167.8 


± 


3.4 


K b 


km s _1 


184.1 


± 


3.4 


CAB 


Rq 


10.99 s 


± 


0.15 


QAB 




0.911 


± 


0.026 


M A 


Me 


3.68" 


± 


0.14 


M b 


M & 


3.35° 


± 


0.14 



a Calculated using the photometric orbital inclination iAB = 88.5 
deg. 

fluxes given bv lClausen et (|l995h . The relative contribu- 
tion to the observed continuum near A5000 was estimated to 
be 0.40, 0.28, 0.16, 0.16 for components A, B, C, and D. Gen- 
erally, when applying a spectral disentangling technique, the 
resulting spectra of the individual components have a very 
high quality, since the S/N-ratio increases roughly as the 
square root of the number of observed spectra involved in 
the calculations (usually several tens), and diminishes pro- 
portionally to the relative flux of the component under con- 
sideration (usually close to 0.5). In the present work a small 
number of spectra are available for a system of 4 spectral 
components. For that reason, the reconstructed spectra, af- 
ter their continuum has been normalized to one, have a S /N- 
ratio lower than the composite observed spectra. In partic- 
ular, for the less luminous stars C and D the effective S/N 
resulted between 50 and 80. 



3 SPECTRAL CHARACTERISTICS AND LINE 
PROFILE VARIABILITY 

The primary component of this subsystem is a typical vari- 
able BpSi star with weak He lines. The spectrum of the com- 
ponent B appears rather normal with strong He and Si lines. 
However, a careful inspection of the high-S/N UVES spectra 
revealed the presence of some asymmetrical and variable line 
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Figure 1. Radial velocity curves for the system AB (upper panel) 
CD (lower panel). Squares, triangles, and small circles are the 
UVES, FEROS, CASLEO observations, respectively. Filled sym- 
bols correspond to the more massive component in each subsys- 
tem (A and C). 

profiles in the spectrum of both stars A and B, which could 
be related to the presence of He, Si, and Mg spots. Since 
a basic assumption for the reconstruction of the component 
spectra is the constancy with phase of the spectral morphol- 
ogy of each star, the disentangled spectra may have small 
artificial features around spectral lines presenting variabil- 
ity. For that reason we preferred, for the analysis of profile 
variability, to use only those spectral lines and orbital phases 
for which the line under consideration is clearly free of con- 
tamination from any other line belonging to the remaining 
components. The first two UVES spectra are particularly 
useful for this purpose since in them both binary systems 
are near opposite quadratures, allowing us to identify clearly 
the four spectral components. 

In Fig. 2 we show the profile of several Si II, Mg II, and 
He I lines, the atomic species that present the most conspic- 
uous lines in the spectra of stars A and B. Both stars A and 
B exhibit line profile variability. Figure 3 shows the spectral 
region around the Mgn doublet at A4481 and the He I line 
A4471. The upper spectra are the two UVES spectra taken 
at quadratures. In the lower spectra, synthetic spectra of C 
and D have been subtracted from the former. The Mgn line 
at A4481 shows the same behavior as silicon lines. In this 
figure we have plotted the observed spectra to demonstrate 
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Figure 2. Line profiles variations for stars A (panel a) and B 
(panel b) near two opposite quadratures: phase 0.34 (left-hand 
part of each panel) and 0.71 (right-hand part of each panel). The 
spectral line profiles shown are, from top to bottom: He I AA4471, 
4026, Sin AA6371, 6347, 5056, 5041, and Mgn A4481. 



that the four components are well resolved, especially star 
B. Consequently, the removal of components C and D does 
not affect significantly the shape of the line profiles. As a 
matter of fact, the profile difference between the two phases 
is evident even in the original composite spectra. In the case 
of star A, the lines are blended with those of component 
C in one of the two UVES spectra, making somewhat less 
confident the analysis of profile variations. 

We note that, even though the reflexion effect due to 
the mutual irradiation of the components would distort the 
line profiles qualitatively in the same way as it is observed 
for Mgn or Sin lines, its expected contribution to the line 
profile variability is much smaller than the observed varia- 
tions. In fact, from the temperature-ratio and the relative 
separation of the components, we estimate that the profile 
variations caused by reflexion effect would be about 5 times 
smaller than the variations observed in Mgn A4481. There- 
fore, the cause of these profile variations is not related to 
proximity effects but to asymmetries in the surface chemi- 
cal distribution. 

Helium lines are weak, especially in star A, so it is not 
possible to detect clear variations in the line profiles from 
the available spectroscopic data. However, we detect that 
the strength of He lines at A4026 and A4471 is variable. 
Other strong He I lines (AA4121, 4922, 5016, 5876) appear 
blended with Fe or Si lines belonging to the companions, 
making them less useful for the analysis of the line profile 
variability. 

For star B the spectral line profiles of Mgn and 
Sill appear strongly variable, as presented in Figs. 2 
and 3. Usually, magnesium does not present large 
equivalent width variations in m agnetic peculiar stars 
l|Leone. Catalano. fc Malarodalll997r i. However, line profile 
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Figure 3. Line profiles variations in components A and B. 
The two upper spectra are the observed UVES spectra for 
MJD 53669.3 and 53670.3, respectively. Labels A, B, C, and D 
mark the position of the Mgn line A4481 for the four compo- 
nents. In the two lower spectra C and D components have been 
removed. The lines Mgn A4481 and He I A4471 are marked for 
components A and B. 



variabilit y of the magnesium do ublet at A4481 has been re- 
ported bv lKuschnig et all (|l999h . who published the first Mg 
map for an Ap star (CUVir). They found that, even with- 
out showing important equivalent width variations, Mgn 
A4481 exhibits notorious profile variations, similar to those 
observed in AO Vel for the components A and B. The pro- 
file shape suggests that Mg is more abundant in the region 
facing the companion star. A more complete study of line 
profiles, that would allow to determine the surface distribu- 
tion of chemical elements, will require a larger number of 
spectra obtained at different phases. 

We can assume that the stellar rotation in the short- 
period eccentric binary AB is synchronized with the angular 
velocity at periastron, so that the rotational period differs 
from the orbital period by a factor 1.163 (P or b = 1.5846 d, 
Prot = 1.3622 d). This fact implies that the stellar surface 
visible at a given orbital phase varies with time, completing 
one cycle in 9.7 days, which corresponds to about 6 orbital 
cycles. A multi-epoch study of AO Vel would be of great 
interest to prove whether the location of stellar spots is fixed 
on the star surface or is related to the position of the stellar 
companion. 

In our spectra of the C and D components the most 
interesting, and rather unexpected fact is the presence of 
spectral lines typical of HgMn stars. The Hgn line A3984 
is present in both stars C and D, while strong Yn and 




M (Msun) 

Figure 4. Mass-radius diagram for the components of the more 
massive binary. The parameters of stars A and B, and their error 
(1-sigma) are represented by the position and size of the elipses. 
Continuous lines are the ZAMS and the isochrone log(age)=7.6. 
Dashed lines, from left to right, are isotherms in the range 
T=12000-14500K, with a step of 500 K. 



Ptn lines are present in the spectrum of component D. Al- 
though comparatively weak, some Mn II lines have been also 
identified in both components. In the reconstructed spec- 
trum of star D we identified Ptn lines at 4046 and 4514 
and more than two dozens Yn lines. No noticeable lines of 
Zrll, Sen, Ban, Pn, Xell, Gall were detected in stars C or 
D. The measured central wavelength of Hgn line at 3984, 
3983.975 A and 3984.08 A for stars C and D, respectively, 
indicates that heavy Hg isotopes are predominant in these 
components. We will return to this point in the next section. 

The weakness of spectral lines for the two less massive 
stars makes impossible to study line profile variations in 
their spectra. 



4 CHEMICAL ABUNDANCES 

The abundances were determined with the synthetic spec- 
trum method usi ng ATLAS9 m odel atmospheres and the 
SYNTHE code jKuruczl Il993h . Stellar temperatures of 
companions C and D were estimated from the correla- 
tion of excitation potentials of Fe lines with abundances 
IjGonzalez. Nesvacil. fc Hubridl2008l) . Microturbulent veloc- 
ities were derived from the correlation of equivalent widths 
of Fe lines with abundances. 

In case of stars A and B, however, it was not possible 
to apply the same procedure since Fe lines are weaker and 
they appear broader because of the higher rotation. Esti- 
mated temperatures for these stars were derived from the ob- 
served masses and radii, i nterpolating in the grid of Geneva 
theoretical stellar models (jLeieune fc Schaereill200ll ) for so- 
lar composition. Fig. 4 shows the position of stars A and B 
in the mass-ratio diagram. Both stars are located close the 
ZAMS. As is shown in this figure, the masses and radii deter- 
mined from light and RV curves correspond to temperatures 
of 13900±500 K and 13200T450K for components A and B, 
respectively. Surface gravity were calculated directly from 
masses and radii. The adopted atmospheric parameters are 
listed in Table 3. 
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Table 3. Atmospheric parameters 



Star 



T eff 

K 



log g 



£ v sin i 

kms" 1 kms 



A 13900 ± 500 4.26 ± 0.03 2 65 

B 13200 ± 450 4.31 ± 0.03 2 60 

C 12000 ± 300 4.3 ± 0.2 2 40 

D 11500 ± 300 4.2 ± 0.2 1 18 



Table 4. Abundances (log[iV '/ 'N (H)]) . Less precise values (error 
= 0.2—0.4) are marked with a colon. 



Element 


Star A 


Star B 


Star C 


Star D 


Sun 


He 


-2.41 


-1.16 


-2.0: 


<-3.0: 


-1.07 


C 


-4.62 


-3.87 


-3.5: 


-3.5: 


-3.48 





-3.71 


-3.56 


-3.7: 


-3.2: 


-3.17 


Mg 


-5.50 


-4.76 


-4.76 


-4.96 


-4.42 


Si 


-3.81 


-4.74 


-4.28 


-4.84 


-4.45 


P 


-5.7: 








-6.55 


s 


-4.67 


-4.67 


-4.67 




-4.67 


Ca 










-5.69 


Ti 


<-7.5 


<-8.0 


-7.0: 


-6.37 


-6.98 


Cr 




<solar 


-6.3: 


-5.87 


-6.33 


Mn 








-5.65 


-6.61 


Fe 


-5.09 


-4.99 


-4.25 


-4.56 


-4.50 


Ni 








<-6.8 


-5.77 


Sr 








-8.27 


-9.03 


Y 








-7.15 


-9.76 


Pt 








-5.94 


-10.2 


Hg 






-5.05 


-5.70 


-10.9 



The results of the abundance analysis are presented in 
Table 4. For comparison we list in the la st column the solar 
abund ances adopted from the work by iGrevesse fe SauvaJ 
(1 19981 ). The abundances of elements not listed in this table 
were assumed to be solar. The complete list of the spec- 
tral lines used for abundaces calculation is included in the 
Appendix A. 

Star A exhibits an overabundance of Si by 0.64 dex and 
underabundance of He by -1.3±0.3dex. This explains why 
He lines are stronger in the star B, even when A is hot- 
ter than B according to the depths of photometric minima. 
Star B appears similar to a normal late B-type star, having 
the abundances of helium and silicon close to solar values. 
However, Fe seems to be slightly underabundant in its atmo- 
sphere. No line of Ti has been observed either in star A or 
B, indicating that Ti is underabundant in both stars. Mg is 
underabundant in star A by —1.1 dex. Similar v alues are fre- 
quen tly observed in magnetic Ap and Bp stars (|Leone et al.l 
1 19971 ). 

We note that the presented abundances in stars A and B 
should be taken with caution, since the detected line profile 
variations suggest a non-uniform chemical distribution of at 
least a few elements. 

The subsystem CD shows abundances which are typi- 
cal for HgMn stars. Star C shows a strong overabundance 
of Hg by 5.8 dex. Other metals like Fe, Ti, Mg, and Si show 
normal solar abundances in star C. Abundances of C, O, 



Ga, Sc, and Sr seem also to be solar. Helium lines are only 
marginally detected in star C and the corresponding abun- 
dance is consequently rather uncertain (about 0.4 dex). It is 
clear, however, that He is underabundant in star C. 

Star D exhibits overabundances of Hg by 5.2 dex, Y by 
2.5 dex, Pt by 5.2 dex, and Sr by 0.8 dex, relative to the solar 
values. Ga and Sc lines are not observed. 

The uncertainties of the abundance determination, typ- 
ically 0.1-0.2, are somewhat larger than in usual analyses 
even though the observed spectra have high S/N-ratio and 
high resolution. Three different reasons contribute to the 
degradation of the accuracy. First, the effective S/N-ratio 
of the normalized spectra of the individual components is 
comparatively low, since the line strength is diluted by the 
continua of the 4 components. In the case of stars C and D, 
for example, the S/N is diminished by a factor 6. Second, 
the disentangling process using a small number of observed 
spectra might generate small fake features, especially close 
to strong variable lines. Finally, the uncertainty of the line 
equivalent widths involves the uncertainty in the relative 
light contribution of each component to the total light of 
the system, i.e. the normalization scale factors adopted from 
the photometric data. We note however that, the chemical 
peculiarities mentioned above are much larger than the un- 
certainties, and the classification of stars C and D as HgMn 
stars is absolutely out of doubt. 

In Fig. 5 we present the synthetic spectrum of star 
D in the spectral region containing Ptll4046, Ptn4061, 
Sr II 4077, and several Fe and Cr lines. As it is visible in this 
figure, Ni is underabundant in this star. Fig. 6 shows the syn- 
thetic spectra for the regions of the Hgll line at A3984. Even 
when the rotational broadening is significant, some informa- 
tion about the isotopic composition can be derived from the 
reconstructed spectra of components C and D. With this 
aim the line profile was fitted with a synthetic spectrum in 
which Hg is assumed to be a mixture of the lightest (196) 
and the heaviest (204) isotopes. As it can be seen in the fig- 
ure, the best fit suggests that the heaviest isotope is by far 
the most abundant in both components of the less massive 
subsystem. 



5 DISCUSSION 

Using high-S/N high-resolution spectra obtained with 
UVES, we have determined chemical abundance for all four 
components of the multiple system AO Vel. Further, line 
profile variations were detected in the spectra of the two 
most massive stars belonging to the AB system. 

Since the multiple system AO Vel is formed by four 
gravitationally bound stars, we can assume that they have 
the same original chemical composition and the same age. 
All components are close to the ZAMS and have already de- 
veloped different chemical anomalies in their atmospheres, 
depending on their temperature and probably due to the 
membership in binary systems. The observed masses and 
radii for the components of the eclipsing pair suggest that 
the system age is less than about 40 million years. 

The connection between HgMn peculiarity and mem- 
bership in binary and multiple systems seems to be sup- 
ported by two observational facts. On the one hand there is 
a high frequency of binaries and multiples among HgMn 
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Figure 5. Comparison of the reconstructed spectrum of the star D and the synthetic spectrum (gray line, red in online version) in the 
spectral region around Pt and Sr lines. In the calculation of this synthetic spectrum, Ni abundance was assumed to be solar so that the 
large Ni underabundance is better appreciated. 
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Figure 6. Observed spectra (thin line) and synthetic spectra 
(heavy line) for stars C and D in the region of the Hgo line 
A3984 and the Yo line A3982. The isotopic composition for Hg 
was assumed to be 196=0.9%, 204=99.1% for starC and 196=4%, 
204=96% for star D. 



stars (|Hubrig. Ageorges. k, Scholleil 120081 ). On the other 
hand, in some binaries the location of chemical spots seems 
to be related with t he positi on of the binary com panion 
i|Hubrig et al.l l2006al . see also ISavanov et al.l 1201391 ). These 
studies indicates that the presence of a companion can play 
a significant role in the development and distribution of 
chemical anomalies on the surface of HgMn stars, proba- 
bly due to proximity effects on the atmospheric temper- 
ature or due to the magnetic field induced in a close bi- 
nary system. The role that magnetic fields possibly play in 
the development of anomalies in HgMn stars, which mostly 



appear in binary systems, has never been critically tested 
by astrophysical dynamos. Recent magnetohydrodynamical 
simulations revealed a distinct structure for the magnetic 
field topology similar to the fractured elemental rings ob- 
served on the surface of HgM n stars (see e.g. Fig. 1 in 
Hubrig, Gonzalez. & The combination of differ- 

ential rotation and a poloidal magnetic fi eld was studied nu - 
merically by the spherical MHD code of iHollerbachl (2000). 
The presented typical patterns of the velocity and the mag- 
netic field on the surface of the star may as well be an indi- 
cation for element redistribution on (or in) the star. 



Generally He and Si variable Bp stars possess large- 
scale organized magnetic fields which in many cases appear 
to occur essentially under the form of a single large dipole 
located close to the centre of the star. The magnetic field is 
usually diagnosed thr ough observat i ons of circular polariza- 
tion in spectral lines. iHubrig et al.l (|2006bl ) used the multi- 
mode instrument FORS 1 mounted on the 8 m Kueyen tele- 
scope of the VLT to measure the mean longitudinal magnetic 
field in the system AO Vel. No detection was achieved with 
the single low-resolution measurement (R=2000) resulting 
in (£? z ) = 80±64G. We should, however, keep in mind that 
the observed spectropolarimetric spectrum is a composite 
spectrum consisting of four spectropolarimetric spectra be- 
longing to A, B, C and D components. Each of the com- 
panions can possess an individual magnetic field of different 
geometry, strength and polarity. Since the inhomogeneous 
elemental distribution of Si and Mg is also observed on the 
surface of the component B, the presence of a large-scale 
organized magnetic field is quite possible in this component 
too. The components C and D show the peculiar ity char- 
acter typical of HgMn stars. IHubrig et al.l (|2006bl) showed 
that longitudinal magnetic fields in HgMn stars are rather 
weak, of the order of a few hundred Gauss and less, and 
the structure of their fields must be sufficiently tangled so 



that it does not produce a strong net observable circular 
polarization signature. 

Thus, the non-detection of a magnetic field in the sys- 
tem AO Vel can be possibly explained by the dilution of the 
magnetic signal due to the superposition of four differently 
polarized spectra. On the other hand, the magnetic field of 
each component should be detectable with high resolution 
spectropolarimeters making use of the Zeeman effect in in- 
dividual metal lines. 

The time scale for the peculiarities to be developed re- 
main unclear since the results of studies of evolutionary sta- 
tus of chemically peculiar stars of different type show some- 
what contradictory results. The early work by lAbtl (|l979h 
suggested that the frequency of peculiar stars of Si and 
HgMn groups increases with age and particularly there ex- 
ist no peculiar star on the ZAMS. In a more recent pho- 
tometric search for peculiar s tars i n five young clusters, 
IPaunzen. Pintado, fc Maitzenl (|2002h also concluded that 
the CP phenomenon needs at least several Myr to start being 
effective. On the other hand, some chemically peculiar stars 
has been detected within very young associations or even 
in star forming regions like Ori OBI dAb t fe Levatoj 19771: 
IWoolf fc Lambertll 19991). Lupus 3 (|Castelli fc HubridboOTf) . 
and L988 dHerbig fe Dahmll2006l) . The results of the present 
paper show the existence of coeval stars with chemical pecu- 
liarities of Si and HgMn types, close to the ZAMS indicat- 
ing that the age threshold for these peculiarities is similar 
for both subgroups of chemically peculiar stars. Our study 
supports the idea that these chemical peculiarities originate 
quite soon after the star formation. 

It is noteworthy, that the subsystem CD, the less mas- 
sive binary of the AO Vel, shows several similarities with 
the eclipsing binary ARAur belonging to the Aur OBI 
association. Both systems are formed by stellar compo- 
nents very close to the ZAMS or even in the pre-MS stage 
ijNordstrom fc Johansenl 1 19941 : iGonzalez et all l2006h . and 
they have almost the same orbital period, 4 d .13 and 4 d .15 
for ARAur and AO Vel CD, respectively. In both systems 
the primary is a HgMn star, but while in the AO Vel CD sys- 
tem the secondary companion shares the same peculiarity, 
in AR Aur the secondary is a normal star (|Khokhlova et all 
1995). In fact, the primary component of ARAur should 
be compared with star D of the system AO Vel, since they 
have a similar effective temperature. This similarity in the 
physical and orbital properties seems to have been trans- 
lated into the development of similar pattern of chemical 
abundances: notorious overabundance of Hg and Pt, over- 
abundance of Sr, Y, and Mn at the level of 1-3 dex, slight 
overabundance of Cr and Ti, and the abundance of Si and 
Mg close to the solar abundance, or slightly subsolar. We se- 
lected a few additional HgMn spectroscopic binaries with ef- 
fective temperatures close to 11,000 K: HD 32964, HD 89822, 
HD 173524, and HD 191110. In Fig. 7 the abundances of var- 
ious elements in AR Aur and selected binaries are compared 
with those of the AO Vel D companion. The abundances are 
presented for HgMn stars member of double-lined spectro- 
scopic binaries with periods between 4 and 12 days and hav- 
ing effective temperatures in the range 10 700-11 700 K. For 
HD 173524 and HD 191110 we plotted the abundances de- 
termined in bot h components. The abundances values hav e 
bee n taken from Adelmanl dl994 ^ Catanzaro fc Letol [|2004ri , 
and lCatanzaro. Leone, fc Letol ( 20031 ). It is remarkable, how 
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Figure 7. Abundances of the AO Vel D companion 
(filled squares) compared with published abundances for an- 
other seven binary components with similar effective tem- 
peratures: HD32964B (circles), HD34364A (open squares), 
HD89822A (pentagons), HD 173524A (stars), HD 173524B (tri- 
angles), HD191110A (plusses), and HD191110B (crosses). The 
dashed line represents the solar abundance. 



much similar the chemical composition of these stars is. We 
note, however, that the abundance distribution on the stel- 
lar surface of HgMn stars is probably inhomogeneous, and 
consequently, the chemical composition derived from spec- 
troscopic observations that do not cover the rotational cy- 
cle, should be taken just as indicative values until accurate 
abundance values obtained using Doppler imaging technique 
become available. The spotted character of magnetic Bp-Ap 
stars is well known, but also in the case of HgMn stars the 
presence of chemical spots or belts on the surface is not 
unco mmon (|Hubrig. Gonzalez, fc Arlt] 120081 : ISavanov et all 
l2009h . 

As already mentioned above, the atmospheric chemical 
composition of ARAur exhibits a non- uniform surface 
distribution with a very in teresting pattern rel ated to the 
position of the companion ([Hubrig et a l. 2006a). A similar 
study in AO Vel would be worthwhile to verify whether 
the elemental distribution on the stellar surfaces of C and 
D components show a similar behavior. Furthermore, an 
intensive spectroscopic campaign that would allow the 
reconstruction of chemical maps for all four components 
should provide important information for the proper 
understanding of the origin of chemically peculiar stars. 



APPENDIX A: LIST OF SPECTRAL LINES 

TablelAl lists all the spectral lines used for abundance deter- 
mination. The lines marked with an asterisk are blends. The 
abundance was determined using the fitting by a synthetic 
spectrum. 

We thank to Charles R. Cowley and Rainer Arlt for 
helpful discussions. 
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Table Al. Spectral lines used for abundace determination. 
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Table Al - continued 
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Table Al - continued 
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YII 3982.592 -0.493 NIST3 -7.45 

4235.727 -1.509 NIST3 -6.60 

4309.620 -0.745 NIST3 -6.95 

4883.684 +0.071 NIST3 -7.25 

4900.120 -0.090 NIST3 -7.25 

5200.406 -0.579 NIST3 -7.45 

5205.722 -0.342 NIST3 -7.15 

PtII 4046.443 -1.190 ENG -5.94 

4061.644 -1.890 ENG -5.94 

4288.371 -1.570 ENG -5.94: 

Hgll 3983.890 -1.510 NIST3 -6.55 

NIST3: NIST atomic spectra Database, version 3 at http://physics.nist.gov 

FW0 6 :lFuhr fc Wiesd ||2006[) 

PTP : IPickering. Thome. &: Perej l|2002h 
LA 
KP 



Lanz fc Artru 



< 1985ft 



TCurucz fcPevtremannl dl975h 
HI : iHibbertl | |l988n 

J07: Johannson (2007), private communication 



K03Cr:Kurucz (2003), http://kurucz.harvard.edu/atoms/2401/gf2401.pos 
K03Mn:Kurucz (2003) ,|http://kurucz.harvard.edu/atoms/2501/gf2501.pos| 
K03Ni:Kurucz (2003),|http://kurucz.harvard.edu/atoms/2801/gf2801.pos| 



K09: Kurucz (2009), http://kurucz.harvard.edu/atoms/2601/gf2601.pos 
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